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Abstract 

Numerical modelling is performed for extreme relativistic parallel 
shocks with upstream Lorentz factor T = 50. Assuming the scattering 
is either large angle or over pitch angles ~ r -1 , spectral flattening and 
shock aaceleration speed-up is found. The energy gain for the first 
shock cycle is ~ T 2 . The likely output from relativistic shocks due to 
the infall from the accretion disk to the AGN black hole is computed. 
Neutrinos from proton-gamma interactions may be detectable with 
planned neutrino telescopes. 



1 Diffusive Shock Acceleration 

Relativistic plasma flow velocities are found in many astrophysical objects such 
as in AGN Jets and their Central Engines where transient shocks are thought 
to be formed and diffusive particle acceleration could potentially occur. The 
gamma flows at these sites are approximatelly 5-10 Lorentz Factors. Gamma 
Ray Burst (GRB) plasma flows also appear to have ultra-relativistic velocities 
where 7 ~ 100 - 1000. Particularly Waxman (1995) and Vietri (1995,1997 and 
1998) have noted (for the case of CRBs) the fact that diffusive acceleration at 
non-relativistic shocks is not fast enough to produce the Highest Energy Cosmic 
Rays (HECR) observed. Vietri (1995) especially, pointed out that simulations 
of Quenby & Lieu (1989) were the first to recognize a considerable enhancement 
in the acceleration rate of a factor~13, for relativistic shocks of 7~3. 

Our current investigation for highly relativistic plasma flows, shows a dra- 
matic acceleration rate enhancement, and since the Lorentz boosting scales as j 2 
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(especially in the first cycle) this effect is important for the Vietri's HECR the- 
oretical predictions. It is known that the chief characteristics of diffusive shock 
acceleration operating in one dimensional astrophysical non-relativistic shocks 
are the 7~ 2 differential spectra for particle Lorentz factor 7 when the shock 
compression ratio is 4 and the analytically derived acceleration time constant 
is given by 



where, A is mean free path, u is flow speed and 1 and 2 refer to upstream 
and downstream. For relativistic flows (eg. AGN jets, T ~ 5 — 10 and GRB fire- 
balls T ~ 100 — 1000) the fractional energy increase for a single shock crossing 
is ~ (1 — V 2 /c 2 )~ - 5 where V is the relativistic difference in flow velocities 
across the shock, but the analytical estimate of cycle time does not change from 
the non-relativistic case and so we expect a speed-up of acceleration relative 
to the non-relativistic estimate. However, because of the large anisotropies 
due to the relativistic flows, a computational approach and careful simulations 
are necessary in the high 7 regime. Several studies have been made by few 
authors, but still, there is a controversy about the shape of the spectrum, the 
rate of energy gain and the time constant for acceleration. However, they are 
crucial to models of AGN Jet and Central Engine cosmic ray acceleration as well 
as the GRB particle production and acceleration. In this paper we present a 
computational study of relativistic diffusive shock acceleration for up to 7~ 50 
and try to investigate the acceleration behaviour at these highly relativistic 
regimes. 

1.1 Numerical Simulations 

In these simulations we study the acceleration of particles by a numerical Monte 
Carlo approach. A guiding centre approximation is used for the particle's propa- 
gation in a ID geometry and the shocks assumed to be "parallel" , either because 
that is the field configuration or turbulence removes "reflection" at the inter- 
face. We keep a compression ratio of 4 in order to be able to be able to compare 
directly with non-relativistic situations and consider large angle isotropic scat- 
tering in the plasma rest frames according to, 
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Figure 1: Energy Spectrum (L), Mean Gamma Gain versus Crossing (R) 

for parallel mfp A, proportional to rigidity and pitch angle 9. Gallant & Achter- 
berg (2000) distinguish a small angle scattering regime where 9 ~ for scat- 
tering models where the field is either uniform or consists of cells of randomly 
orientated uniform field. For the time being we are not in this limit, but con- 
sider the cases where the particle coming upstream from the shock moves a 
distance ~ A with only a small net total deflection, during which time it must 
be travelling close to the shock normal, otherwise it could not penetrate up- 
stream against the relativistic convective flow. Then a turbulent field or flow 
discontinuity is encountered and deflection takes place in a time ~ r 9 /c where 
r g is the gyroradius. Provided A > r g T, significant deflection may take place 
before the particle is swept back to the shock. Lorentz transforms are used be- 
tween frames on shock passage and the results recorded just downstream in the 
shock frame. In this numerical investigation and in order to keep the statistics 
of the results as reliable as possible, a particle 'splitting' technique is used. The 
particles injection is upstream, where approximately 10 6 particles have been 
used for the completion of the code runnings. 

Results for isotropic scattering with an upstream T ~ 50 are shown in figure 
1(L) where the energy scale is measured by the particle Lorentz factor, together 
with the best fitting power law. Clearly other factors, especially trapping, limit 
the actual high energy extent of the spectrum. The plateau-like structure corre- 
sponds to successive cycles of acceleration. Our calculations show that roughly 
80% of all particles are lost downstream each cycle, because of the high plasma 
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Figure 2: Acceleration Speed-Up, Isotropic Scatter (L), Pitch-Scatter (R) 



flows. The mean energy gain per cycle on cycle 1 for isotropic scattering is 80 
and on cycle 3 is 3200, to be compared with the square of the upstream flow 
gamma factor r^ p =2500. This is seen in figure 1(R) where injection is at 7 = 50. 
Figure 2(L) shows the ratio of the computational to the non-relativistic, analyt- 
ical acceleration time constants as a function of T up , for isotropic scattering and 
figure 2(R) for small-angle scattering. For small- angle scattering, the mean free 
path to scatter ~ n/2 is used in the analytical expression. As it may be seen, a 
substantial 'speed-up' is noticed in both cases and along with our results of the 
very large energy boosting in consecutive shock cycles, which establishes the 
allowed upper energy limit in GRB fireballs, Vietri's theory is been supported 
regarding the predictions of the origin of Ultra High Energy Cosmic Rays from 
GRB fireball sites. 



2 AGN Accretion Shock Model 

A preliminary application of the speed up and spectral hardening confirmed 
above has already been given by Battersby et al (1995). Central accretion shocks 
for infall to AGN black holes must be sporadic because stable shocks occur at 
a given position for a very specific plasma angular momentum and internal en- 
ergy (Molteni, Lanzafame & Chakrabarti 1994). We expect quasi-perpendicular 
shocks to form out to 100 Schwarzchild radii because quasi-keplerian plasma mo- 
tion should be taken into account but is neglected in Battersby et al (1995). We 
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estimate the lifetimes as given by free-fall and get 'de Hoffmann- Teller' frame 
speeds > 0.5c Using the above calculations of speed-up and spectral slope and 
allowing a realistic increase in the mean free path above cyclotron radius yields 
the acceleration time and spectrum. The input particle total energy and mag- 
netic field is estimated assuming Lx- ra y — 0.05L^^ is in equipartition with 
magnetic energy and particle energy, B oc R^ 1 and the Stecker-Salamon photon 
spectrum is used. Proton scattering in the field and convection into black hole 
is followed for radial infall and quasi-keplerian plasma motion. The secondary 
cascade following initial acceleration and propagation is followed. High energy 
protons cool preferentially on photons while low energy protons cool preferen- 
tially on background protons, pion decay into 7, v and e ± components and 
also a 7 cascade, synchrotron photon production and 7 — 7 interactions result. 
Battersby et al (1995) concentrate on neutrino production from p gamma in- 
teractions. These interactions are more important than pp at high energies. 
They find that a 1 km 2 high energy neutrino telescope is really required to have 
a good chance of detecting fluxes from a source like 3C273 or the background 
from a large number of sources. 

3 Conclusions 

A Monte Carlo numerical investigation has been reported on the test particle 
limit of parallel diffusive shock acceleration. Very high gamma flow astrophysi- 
cal plasmas have been used up to 7„ p ~ 50 which are relevant to AGN Jet and 
Central Engines shocks as well as the GRB ultra-relativistic enviromments. A 
dramatic 7 2 (especially for the first cycle) energy boosting has been observed 
and a significant acceleration rate efficiency has been reported. Especially, both 
these results support (and been supported by) Vietri's and Waxman's models 
and theoretical predictions on the problem of HECR origin from GRBs. This 
crucial issue though, of ultra-relativistic diffusive shock acceleration, needs fur- 
ther investigation as many other parameters are needed to be included within 
the simulation codes in order for the model and the results to be as realistic as 
possible in connection with the observations and the theoretical predictions. 

Much of the secondary output from the AGN originates from sporadic shocks 
near 100 Schwarzschild radii, pp interactions are less important than p7 in neu- 
trino production except at low energy. 7s from synchrotron dominate 7s from 
7rs up to 1 TeV, but the GeV to TeV 7 flux is a factor 100 less than observation 
of 3C273 output, suggesting some "Beaming" or amplification mechanism is 
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required if central accretion shocks are the main source. We expect significant 
fluxes from < 100Rs c h if the short period fluctuations observed are related to 
black-hole physics. The 100Rs c h proton output < 0.1 primary proton output 
> 1 PeV from the central shocks. The neutron output < 0.01 proton primary 
output. 
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